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Feasibility of approximating spatial and local entanglement in long-range interacting 

systems using the extended Hubbard model 
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We investigate the extended Hubbard model as an approximation to the local and spatial entan- 
glement of a one-dimensional chain of nanostructures where the particles interact via a long range 
interaction represented by a 'soft' Coulomb potential. In the process we design a protocol to cal- 
culate the particle-particle spatial entanglement for the Hubbard model and show that, in striking 
contrast with the loss of spatial degrees of freedom, the predictions are reasonably accurate. We 
also compare results for the local entanglement with previous results found using a contact inter- 
action 1] and show that while the extended Hubbard model recovers a better agreement with the 
entanglement of a long-range interacting system, there remain realistic parameter regions where it 
fails to predict the quantitative and qualitative behaviour of the entanglement in the nanostructure 
system. 
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I. INTRODUCTION 

The properties of quantum dot-based nanostructures 
make them possible candidates as hardware for the ma- 
nipulation of quantum information 2!-f6l| . Chains of quan- 
tum dots (QDs) have been proposed to transfer quan- 
tum information Q, and to generate, distribute, and 
freeze entanglement @, Q • Entanglement is in fact con- 
sidered a fundamental resource for quantum information 
processes and is increasingly gaining the attention of the 
condensed-matter community as a probe for delicate phe- 
nomena, such as quantum phase transitions |10 |. 

The Hubbard model '1V\ is a simplified model of itin- 
erant interacting fermions with positions discretized to 
a lattice, and usually only includes interaction within a 
lattice site, while extended Hubbard models may include 
long-range interactions. The calculation of the entan- 
glement in a realistic system of many fermions, such as 
electrons trapped in a nanostructure, is usually compu- 
tationally too demanding. In this respect, proving the 
Hubbard model - or a variant of the Hubbard model - 
accurate as an approximation to the entanglement of a 
realistic many-fermion system would open the possibility 
of using density-functional theory techniques to calculate 
the entanglement of these complex systems: a powerful 
local-density approximation approach is in fact available 
for calculating the entanglement of the Hubbard model 

In previous work 1] the accuracy of the ID Hubbard 
model (HM) as an approximation to the average local en- 
tanglement of two-electrons trapped in a chain of quan- 
tum dots (QDs) was considered. In [y we considered a 



contact interaction between particles and found the HM 
to be accurate for single-site entanglement calculations, 
but the spatial entanglement could not be estimated as 
there was no scheme for such entanglement measurement 
in the HM. 

In this paper we consider a similar system but with 
a more realistic long range particle-particle interaction, 
which we model as a 'soft' Coulomb potential. The accu- 
racy of the local entan glem ent of the ID extended Hub- 
bard model (EHM) pl.Tli| and the EHM with correlated 
hopping (EHM+CH) as an approximation to that of the 
QD structure is then appraised. In addition we propose a 
method to calculate the spatial entanglement of the HM 
and EHM and compare this somewhat severe approxi- 
mation with the spatial entanglement of the QD system. 
Surprisingly our results show that even if in the HM and 
EHM cases the number of spatial degrees of freedom are 
reduced to just a few lattice sites, the HM and EHM spa- 
tial entanglement still capture most of the features and 
trends of the spatial entanglement of the quantum dot 
structure, and it does so for both attractive and repul- 
sive interactions. 



II. THE QD-CHAIN TWO-PARTICLE SYSTEM 

We consider a system of two fermions trapped within 
a one-dimensional QD-chain. In effective atomic units, 
the Hamiltonian is 






^QD= E [-,^^+vi^s))+Cuf{K-x,\). (1) 

s—a.b 
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Here v{xs) is the potential used to represent an array of 
regularly spaced, identical square wells, each well repre- 
senting a QD. The chain is symmetric around the origin, 
and defined by the quantities: M the number of wells, d 
the barrier width between two consecutive wells, and w 



and vq the width and depth of each well respectively. We 
define and vary the interaction strength Cu to facilitate 
comparison of the QD system with the HM. The inter- 
action type we consider is either a contact interaction, 
f{\xa — Xb\) = S{xa — Xh}, or a long range interaction of 
the form (see e.g. refs. p^ll6|) 



f{\Xa -Xt\) = 



1 



^iP + {Xa-Xbyy 



(2) 



This is often referred to as a 'soft' Coulomb potential. 
Here we use I = 1 ao, ao the effective Bohr radius. In 
this work we focus on a system of four wells as previ- 
ous results |lj| showed that this was the smallest number 
for which the average local entanglement exhibits a non- 
trivial dependence on Cjj. We calculate the solution of 
the time-independent Schrodinger equation correspond- 
ing to eq. (H)) by using 'exact' diagonalization with a basis 
formed by the eigenfunctions of the non-interacting sys- 
tem. We note that the system ground-state is a singlet 
due to the choice of zero magnetisation. 



III. AVERAGE SINGLE-SITE (OR LOCAL) 

ENTANGLEMENT 

Here we consider the average single-site (or local) en- 
tanglement of the system ground state. This type of 
entanglement is relevant for systems of indistinguishable 
fermions [13, [l3 ■ To this aim we divide our QD system 
into contiguous 'sites', each site centred around a single 
well. The entanglement entropy S of the system is given 
by 



1 ^^ 



(3) 



with Si = —Trprcd,i log2 Prod.i the j-site von Neu- 
mann entropy of the reduced density matrix prod.i- 
By dividing the system into sites and moving to a 
site-occupation basis the reduced density matrix be- 
comes the 4x4 diagonal matrix [13, Il3 Prcd.i = 
diag[P,{n),P^{f},P^{i),P^{0)] ,with P,{a) the probabil- 
ity of double (a =ti), single {a =t or |), or zero {a = 0) 
electronic occupation at site i. 

We calculate the ground-state wavefunction, for an 
even number M of wells, and from that obtain the occu- 
pation probabilities, as described in detail in ref. [jy]. 



IV. COMPARISON WITH HUBBARD MODEL 
VARIANTS 

We consider the Hamiltonian 

+U'^ni^^ni^l + U' '^ fii^crni+i^a', (4) 



where U and U' are respectively the on-site and inter- 
sites interaction strength, t is the hopping parameter and 
i' is the correlated hopping term, c} ^ (ci,cr) are the cre- 
ation (annihilation) operators for a fermionic particle of 
spin a at site i, and fii^cr = cl ^Ci^a is the particle number 
operator. To solve eq. @ we use exact diagonalization in 
the single-site occupation basis {|ti) , It) , li) , |0)} with 
open boundary conditions and an average particle den- 
sity oin = ni+n-f — 2/M . Here n^ is the average density 
of the cr-spin component. Again we calculate the average 
single-site entanglement according to eq. ^ 19l-l2l|. We 
note that for t' — 0, eq.(|3]) describes the EHM while for 
f/' = ? = , the HM is recovered. 

In order to calculate the equivalent of t for the QD 
system, Iqe), we use the same procedure employed in [l| 
for the contact interaction. The corresponding on-site, 
inter-site and correlated hopping for the QD model are 
calculated from 

J Cu f (j>h{Xa)<i>%{.Xb)<i)j{Xa)4'k{xb) J , ,^. 

hzjk = -rr- / j===^=^== dxadxb, (5) 

^ J \/{l^ + [Xa-Xby) 

with Uqd — Illll, U'qjj = Ilrlr, tgu — Illlr. and 
4>L(R) the single-particle ground state of the finite single 
square well potential, but positioned in the left (0l) or 
right {(J)r) well. 

Usually the hopping parameter t is used to rescale C/, 
U' and t', giving the dimensionless interactions U = f//t, 
U' = U' /t and t' = I' jt. For the QD system we obtain 
U = UQo/tQD = 17865C[/, U' = 0.28^7 and t' = 2 x 
lO^^C/ for d = w = 2 ao, and uq — 10 Hartree. 

We first consider the HM. Fig. [l] shows that the HM is 
not as good an approximation when long range interac- 
tions are used. In the QD system with d — 2ao, the 'long 
range' average single-site entanglement for < C/ < 10 
has some of its behaviour captured by the HM but the 
long range interaction induces a higher maximum which 
occurs at a different position, and there is no flex point 
in the curve. The calculation of Uqd and igo using 
the single-particle single square well ground-state may be 
expected to be a more severe approximation for the long 
range than for the contact interaction case. This does not 
appear to be the main reason for the difference shown in 
fig. [1] though, as a scaling of U by fitting UQD/{tQDCu) 
cannot rectify this difference. This suggests that there 
are significant contributions to the interaction beyond 
the on-site repulsion. 

We also investigate the long range interaction case 
when the wells are further apart {d — 4ao): this would 
be expected to reduce the electron density in the barrier 
region, and make the sites better defined and hence the 
system more similar to the HM. We note that our method 
for calculating Iqd becomes in this case too susceptible 
to the noise in the tail of the numerical wavefunction, as 
there is an almost negligible overlap between the single- 
particle wavefunctions centred in the right and left well. 
So this method fails to give a reliable tgn value. Hence 
for d = 4 oo we use tqD as a parameter to generate a 
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FIG. 1; Average single-site entanglement vs interaction 
strength [/ for M = 4. HM (dotted line) and QD system with: 
contact interaction, U = 12344CJ7 and d = w — 2 ao (solid 
line); long-range interaction, U — 17865C(7 and d = w = 2 ao 
(thick dashed line); long-range interaction, U = 8.5 x 10^ Cu, 
d — 4 ao and w = 2 ao (thin dashed line) . For all QD systems 
vo ~ 10 Hartree. 



QD entanglement curve closer to the one representing 
the HM system. With this procedure, the average single- 
site entanglement's maximum is closer to that of the HM 
than for d = 2 ao, but the overall shape of the entan- 
glement curve resembles the d = 2 ao case, as shown in 
fig. [TJ It therefore strongly suggests that much of the 
dissimilarity with the HM is caused by the long range in- 
teraction rather than the proximity of the wells. In fact, 
unlike a contact interaction, a long range interaction does 
not require the particles' density to overlap. To test this 
hypothesis we consider next the EHM and EHM-I-CH 
in which interactions between particles on neighbouring 
sites are included. 

For d ~ 2 ao the entanglement results for the EHM and 
EHM-f CH are indistinguishable on the scale of the plots 
(see fig. [2]). We see in the upper panel of fig. [5] that the 
EHM reproduces better the behaviour of the entangle- 
ment when long range interactions are considered, with 
the most appropriate U' seemingly residing somewhere 
between 0.15C/ and 0.2U. This difference with our es- 
timated value of U' — Q.28U could arise from the use 
of the non-interacting single square well solutions in our 
calculation of U' . 

For attractive interactions, the lower panel of fig. [2] 
shows that the EHM reproduces the single-site entan- 
glement of the long-range QD system fairly well at all 
U values, as well as the HM reproduces the single-site 
entanglement of the QD with contact interaction. For 
repulsive interactions, although the contact interaction 
and the HM embody some of the features of the long- 
range interaction system for the parameters chosen, the 
average single-site entanglement arising from long range 
interactions is lower for U ^ 8, significantly so for very 
large U. This is due to the long range repulsion forcing 
nearly all of the particle density into the outer wells: the 
particle density becomes significantly different from zero 
only in the outer wells where though only single occu- 
pation remains non-negligible. For this case the entan- 
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FIG. 2: Average single-site entanglement vs interaction 
strength U for M = 4. Upper panel: long range interaction, 
d = w ^ 2 ao, U = 17865C(7 (thick dashed line); EHM and 
EHM-HCH (circles) with: U' = 0.15(7 (solid line); U' = 0.2U 
(dotted line); U' = 0.28U (thin dashed line). Lower panel: 
contact interaction, d = w = 2ao,U = 12344C(7 (solid 
line); HM (dotted line); long-range interaction, d — w — 2 
ao, U = 17865Cu (thick dashed line); EHM (thin dashed 
line) and EHM-f CH (circles) with U' = 0.28(7. For aU QD 
systems no = 10 Hartree. 



glement would be expected to be bounded from below 
by 2/M. Fig. [2] shows that in the long range case the 
entanglement does indeed come close to this limit. This 
density configuration is not accessible with the contact- 
type interaction: this interaction could affect the density 
profile to some extent, but could not make the outer wells 
more favourable than the inner wells (see [1[ and fig. [5]) . 
The EHM reproduces fairly well the effect on the den- 
sity of long-range interactions for U < 10, and we see 
that the related entanglement saturates at a value lower 
than the HM for higher values of U. However this satu- 
ration value is still much higher than the entanglement 
values reached by the QD system with long range inter- 
actions, which do not yet saturate even for U as large as 
150. On the other hand, the greater effect on the density 
distribution means that the QD system with long range 
interaction comes closer to achieving the theoretical max- 
imum entanglement for M = 4, S'^q^. — 1.623 (see table 
I in ref. [1[): for long-range interactions Smax = 1-565 
compared with S'max — 1-550 for the contact interaction 
(see fig-[l|)- 

We see in fig- [3] that in the limiting case of d = 0, for 
which our QD model reduces to a single dot of width 
AIw, the EHM is fairly accurate for the entanglement of 
the QD system when \U\ is small- However the maximum 



and U ^ values are less accurate than with d = 2 gq: 
compare the upper panels of figs. [5] and [31 We note that 
for d = the HM and EHM are also substantially less 
accurate for [/ < than in the d = 2 aa case (compare 
lower panels of figs. [2] and [S]). Including the correlated 
hopping term modifies only slightly the entanglement of 
the EHM for C/ > 0, but has a noticeable effect for [/ < 
where it produces a better approximation to the long 
range QD system. 



rable or larger than w. However they all appear to be 
a fair approximation to the average single-site entangle- 
ment in the limiting case of d = and small U suggesting 
that the HM, EHM and EHM-f CH are useful in mod- 
elling systems of dots when d is much smaller than w as 
such systems correspond to small U. Interestingly, this 
may include the case of a single dot when modelled as a 
(small) set of finite partitions (system GaAs'^'^^' in table 
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FIG. 3: Average single-site entanglement vs interaction 
strength U for M = 4. Upper panel: HM (dotted line); long 
range interaction, d = 0, w — 2 ao, U = l.SACu (thick dashed 
line); EHM with U' = 0.54(7 (thin dashed hne); EHM+CH 
with [/' = 0.54C/, t' = 0.047[/ (circles). For all QD systems 
vo = 10 Hartree. Lower panel: same as upper panel but for a 
greater range of interaction strengths. 

The Hubbard model with U < 8 is used to describe 
strongly correlated solids [22], while larger U (up to 
~ 150) can represent systems of optically trapped atoms 
[23| . To estimate the typical values of C/ in a quantum 
dot system, we may use effective atomic units and the 
observation that a system with Cjj = K, well depth vq, 
well width w, barrier width d and minimum interaction 
length I is equivalent to a system with Cu — 1 , well depth 
Vq/ K"^ , well width Kw, barrier width Kd and minimum 
interaction length Kl. Under these conditions, given a 
specific QD system and the well parameters, the on-site U 
can be estimated, as shown in table H] We then note that 
some typical QD parameters such as in !24!] and f23| cor- 
respond to values of U for which the EHM and EHM+CH 
are unable to reproduce the QD entanglement. 

The results in table U suggest that for a long range 
interaction the HM, EHM and EHM+CH are not at all 
suitable to accurately model the average single-site en- 
tanglement of realisable quantum dots when d is compa- 



V. PARTICLE-PARTICLE SPATIAL 

ENTANGLEMENT 

We now investigate spatial entanglement y^, i.e. the 
particle-particle entanglement related to the spatial de- 
grees of freedom. For the QD system the related reduced 
density matrix is infinite dimensional and given by 



/5rcd,spatial(a;,x') = / ^qq (x, Xr )*QD (a:', Xr )dXr (6) 

where ^qd is the QD system wavefunction. 
The von Neumann entropy is then 5'spatiai = 
-Trpt.od,spatiailog2Prod,spatiai- Wc cousidcr the QD 
system with both contact and long-range interactions, 
and compare it with the HM and the EHM (with and 
without the CH term) with M sites. To calculate the 
HM (EHM) spatial entanglement we transform the HM 
(EHM) wavefunction from the site occupation basis to 
the particle basis to give ^hm (ehm)- We move from 
having a superposition of the tensor products of M 
sites, where each site can be in one of four states, to a 
superposition of the tensor products of two particles, 
each of which can be in one of M sites. This mapping, 
for the two-electron system under consideration, is given 
by 



|0)i. 

5^[lta)|;6) 



in),...|o),, 

\ia)\U)l 



\Xi.b) 



(7) 



for states' components with doubly occupied sites, and 
by 



\Xj,a) \Xifi) \ia) Wb) 



(8) 



for states' components with single-occupied sites, where 
antisymmetry in particle exchange has been imposed. 
Here {xi^s}, with i = 1,2...M, and s — a, 6, repre- 
sent the discrete M possible coordinates (site positions) 
for particle s. The spin part of the resulting ground- 
state wavefunction factorises into a singlet, so its entan- 
glement is constant and, as such, not of interest here. 
The spatial part of the wavefunction vI^hm (ehm) can now 
be used to calculate the M x M reduced density matrix 

p"Jd}p"i^P = TrA |*HM (EHM)) (*hm (EHM) | where we 



TABLE I: Estimate of the on-site interaction strength U for different QD systems: GaAs-type systems with reduced mass 
nie// = 0.067me and dielectric constant e = 10.9 (GaAs'^' to GaAs''^'); GaAs-based system with m^ff — 0.067me and e = 12.1 
(GaAs^*' [2JI)j and CdSe-based system with m^ff = 0.45771^ and e = 9.1 25J. For the d = hmiting case, the system is 
physically a single dot of width 44nm but modelled as four partitions of width w = 11 nm. 
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GaAs'i' 
GaAs(2) 
GaAs'^' 
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CdSe [25] 



QD Parameters 

iio(eV) «;(nm) d(nm) I (nm) 



1.5 5.5 5.5 
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model parameters 

~ U 
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trace out the subsystem A corresponding to one of the 
particles. The reduced density matrix is then used in the 
evaluation of the von Neumann entropy ^spatial- 

Due to the lattice discretization, as an approximation 
to vI/qd, ^l/HM (EHM) retains very few of the spatial de- 
grees of freedom, as few as 4 in the system at hand. Yet, 
our results in fig.|4]show that the agreement between the 
EHM spatial and the QD system spatial entanglement is 
surprisingly good, especially for attractive interactions. 
We also find that the HM spatial entanglement is an ex- 
cellent approximation to the QD system with a contact 
interaction for both attractive and repulsive interactions 
(fig. |4]). Interestingly, the QD system with long range 
interaction has a lower entanglement for large positive U 
than the QD system with contact interaction, even if the 
former has stronger. Coulomb-dependent, correlations. 
That occurs because the long range interaction shapes 
the particle density substantially more, which means that 
for large repulsion almost all of the density resides in 
the outer wells, see fig. [51 This results in a triplet-type 
state 16] for which the spatial entanglement is equivalent 
to the entanglement of a maximally entangled two-qubit 
state, i.e., unity. Our data confirm this picture. Interest- 
ingly in this regime the EHM is actually a poorer approx- 
imation than the HM to the long range QD system. The 
EHM has increased Coulomb correlations in respect to 
the HM and hence displays a higher entanglement. How- 
ever the EHM fails to effectively exclude the particles 
from the inner wells and to reproduce the triplet-type 
state (with its lower entanglement) which characterises 
the QD system when U » 10. In fact at U = 100 
the single-particle occupation at the outer (inner) wells 
is 0.96 (0.045) for the QD system with long-range interac- 
tion, while for the EHM the probability of single-particle 
occupation at the outer sites is 0.74 and at the inner sites 
it is still as high as 0.26. This implies that the rank of 
the EHM density matrix does not reduce to the one of a 
two-qubit density matrix. 



VI. CONCLUSION 

We studied the Hubbard, extended Hubbard, and ex- 
tended Hubbard with correlated hopping models as an 
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FIG. 4: Spatial entanglement vs interaction strength U: con- 
tact interaction, d = w — 2 ao, U — 12344C(7 (solid line); 
HM (dotted line); long-range interaction, d = w = 2 ao, 
U = 17865C!7 (thick dashed line); EHM (thin dashed hue), 
EHM-hCH (circles) with U' = 0.28(7. vo = 10 Hartree for all 
QD systems. 




FIG. 5: Density profile for the QD system ai U = 100 with 
A4 = 4, d = w = 2 ao, Vo = 10 Hartree and a contact (solid 
line) or long range (thick dashed line) interaction. 



approximation to the local and spatial entanglement 
of a one-dimensional chain of nanostructures, in which 
trapped particles are interacting via long-range interac- 
tions. We focused on a system comprised of 4 quan- 
tum dots, though the typical trends we found apply to 
longer chains (not shown), as the long-range interaction 
will eventually force the density into the outer wells. 

For well-separated QDs the EHM reproduces fairly well 
the average-single site entanglement for attractive long- 



range interactions and < U < 10, though with a max- 
imum positioned at a weaker interaction strength. In 
contrast with the contact interaction case [l|, for fong 
range interactions and U >10 the EHM completely fails 
to reproduce the single-site entanglement, with an error 
which rapidly increases up to ^ 150% at U = 500 (not 
shown). We emphasise that [/-values as high as 5700 
do indeed correspond to some of the typical parameter 
ranges for QDs, so care must be taken when trying to 
model nanostructure systems with the HM and EHM. 

We note that ref. [23] found the HM not to be an 
accurate approximation when considering the exchange 
coupling in two-electron double quantum-dots modelled 
with a parabolic confining potential and a Coulomb in- 
teraction in two dimensions or an effective model in one 
dimension. However the HM did at least qualitatively 
reproduce the behaviour of the exchange coupling as the 
intcr-dot distance was varied, while, interestingly, the 
EHM did not. In contrast we found the EHM to be a 
better approximation than the HM to the average single- 
site entanglement of the QD-chain system, though with 
the limitations described above. 

We considered the limiting case of interdot barrier 
d = 0. In this case our model describes a single QD, 
which is divided into an arbitrary number M of parti- 
tions. We then considered the average entanglement of 
one of these partitions with the others. We found that 
this corresponds to small U for physically realisable dots 
and in this region the HM and EHM accuracy as an ap- 
proximation to the entanglement were fairly good. 

We also studied the case of attractive interaction: for 
well-separated QDs the average single-site entanglement 
was well approximated by both the HM and the EHM for 
any interaction strength. However as the interdot barrier 
width decreases, the HM and EHM fail to reproduce the 
entanglement of the QD system at medium and large 
interaction strength, quickly reaching a discrepancy of 
~ 40% for the limiting case d = 0. 

We designed a procedure to calculate the spatial en- 
tanglement from the HM and found that the spatial en- 
tanglement associated with the HM was a very good ap- 
proximation to the spatial entanglement of the four QD 
chain when particle-particle contact interaction was con- 



sidered. This was surprising as we lose most of the spatial 
degrees of freedom when using the HM to calculate the 
spatial entanglement: we consider one position per site 
only - four positions in total for the specific case anal- 
ysed in this paper - compared with an infinite number 
for any QD system. Yet our results show that the perti- 
nent information seems to be retained. The QD system 
with long range, repulsive interaction was less well ap- 
proximated by the HM, as could be expected, but inter- 
estingly the EHM performed worse still. However their 
errors at large U, ~ 13% and ^ 33% respectively, were 
considerably smaller than for the average single-site en- 
tanglement of the long-range QD system. We attributed 
the worse performance of the EHM to the long range in- 
teraction causing most of the electron density to reside 
in the two outer wells so essentially giving the maximum 
entanglement of a two-qubit system. The EHM fails to 
modify the density to this extent, and hence to lower 
the entanglement when U >> 10. We find that for the 
analysed QD system the correlated hopping term is very 
small and produces appreciable effects on the entangle- 
ment only in the limiting case of d = with attractive 
interactions. 

Our results show that when dealing with nanostruc- 
tures care must be taken in assessing the validity of the 
HM or EHM as an approximation to the system be- 
haviour, as for some of the experimentally relevant pa- 
rameter range these approximations fail badly to repro- 
duce the local entanglement. Surprisingly, our calcula- 
tions show that for QD chains the essential properties of 
the spatial entanglement are retained when discretizing 
the spatial degrees of freedom to a handful of points, and 
in this case the HM and EHM remain fair approximations 
at all parameter values. 
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